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ABSTRACT 

A  theoretical  study  of  the  vibrational  relaxation  of  the  adbond  between  a 
physlsorbed  molecule  and  a  crystalline  substrate  is  presented.  The  bond 
between  the  admolecule  and  the  substrate  is  described  as  a  one -dimensional 
Morse-oscillator.  Using  standard  perturbation  theory,  expressions  for  the 
relaxation  of  the  vibrational  adbond  due  to  its  interaction  with  the  substrate 
lattice  vibrations  are  derived.  Both  energy  relaxation  and  pure  dephasing  are 
considered.  A  numerical  evaluation  of  the  resulting  expressions  for  the 
relaxation  constants  is  given  for  a  Debye  spectrum  for  the  phonons  and  for  a 
surface  phonon  spectrum.  The  relative  importance  of  pure  dephasing  and  energy 
relaxation  in  determining  the  absorption  llnewldth  is  derived  as  a  function  of 
the  fundamental  frequency  of  the  active  mode. 
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1.  INTRODUCTION 


The  effects  of  vibrational  excitation  of  adsorbed  molecules  by  an  IR 

laser  has  been  the  subjeet  bf.many  experimental  and  theoretleal  studies  over 
«  ' 

the  past  decade.  Primarily,  the  Interest  in  laser  excitation  of  an  adspecles 

2 

arose  because  of  the  prospect  of  influencing  surface  chesical  reactions  and 

the  prospect  of  selective  laser-induced  desorption  without  heating  (l.e., 

damaging)  the  substrate.^  Later  it  was  shown  that  substrate  heating  via  the 

adsorbate  of  an  otherwise  transparent  crystal  sight  quite  well  be  the  most 

D-A 

important  process  occurring. 

Besides  using  the  laser  to  Influence  particular  processes,  it  can  also  be 
used  as  an  analytical  tool.  Best  known  is  the  measurement  of  the  absorption 
line  profile  and  the  line  width.  More  modern  examples  are  echo  experiments 
and  time-resolved  transmission  spectroscopy.  With  all  these  techniques  some 
information  can  be  obtained  about  the  adbond  or  adsorbate  and  Its  Interaction 
with  the  substrate. 

Surface  reactions  are  mainly  dependent  on  the  extent  to  which  a 
particular  vibrational  mode  can  be  excited  and  the  speed  with  *ilch  relaxation 
ooours.  Processes  in  which  the  populations  of  the  vibrational  levels  are 
ohanged,  due  to  an  energy  exchange  with  the  reservoir,  are  generally  called 
energy  relaxation  or  T^ -processes. 

The  absorption  llneshape  is  determined  by  the  Fourier  transform  of  the 
autocorrelation  function  of  the  dipole  operator  of  the  adbond.  In  particular, 
it  Involves  the  time  evolution  of  the  off-diagonal  elements  of  the  reduced 
density  operator.  In  the  most  simplest  approximation  the  absorption  line  is  a 


Lorentzian,  with  a  width  determined  by  two  processes.  The  first, 
inhomogeneous  broadening,  arises  due  to  inhooogeneities  in  the  environment  of 
the  active  molecules.  For  adsorbates,  it  can  have  contributions  from 
differences  in  adsorption  sites,  roughness  of  the  surface  and  other  causes. 
Because  the  magnitude  of  these  irregularities  la  not  well  known,  inhomogeneous 
broadening  is  difficult  to  analyze  in  a  first  principles  way.  A 

o 

phenomenological  approach  was  recently  given  by  Gortel  et  al. 

The  second  mechanism,  homogeneous  broadening,  again  consists  of  two 
parts:  energy  relaxation  (or  -  processes)  and  pure  dephasing  (or  T|- 
proc esses) .  The  latter  is  the  collection  of  all  relaxation  terms  which 
contribute  to  the  time  evolution  of  the  off-diagonal  elements  of  the  density 
operator,  but  tdilch  do  not  change  the  populations.  In  the  simplest 
approximation  the  total  homogeneous  width  is  given  by. 


1  2.1 

T  "  T  T* 
2  1  2 


(1.1) 


Each  of  these  three  relaxation  times  can  in  principle  be  determined 
experimentally.  For  adsorbates,  the  energy  relaxation  tine  was  determined  for 
chemisorbed  OH  on  allies  by  tine-resolved  transmission  spectroscopy. 10,11  The 
same  experiments  for  the  lower-frequency  adbond  node  of  physlsorbed  systems 
have  not  been  done  yet.  One  problem  there  is  the  lack  of  an  appropriate  laser 
for  these  frequencies.  Echo  experiments  give  the  total  dephasing  tine 
(Eq.  (1.1)).  Together  with  the  previous  mentioned  time  resolved  spectroscopy 
measurements,  they  can  be  used  to  determine  the  pure  dephaslng.  Finally, 


simple  measurement  of  the  absorption  llnewldth  gives  the  eoablnatlon  of 
homogeneous  and  inhomogeneous  broadening.  Coabined  with  the  result  for  T1  and 
T2«  it  gives  the  contribution  of  inhomogeneous  broadening. 

The  question  we  want  to  address  in  this  paper  is  what  the  relative 
iaportance  is  of  energy  relaxation  and  pure  dephasing,  respectively,  to  the 
hooogeneous  llnewldth. 

Experimentally,  we  know  from  the  -measurements  by  Heilweil  et  al. 10,11 

12 

and  previous  determinations  of  the  total  llnewldth  that  the  latter  is  not 
due  to  Tj -processes  only.  Whether  the  discrepancy  has  to  be  attributed  to 
pure  dephaslng  or  to  Inhomogeneous  broadening,  cannot  be  deduced  from  these 
experiments.  Measurements  of  the  total  llnewldth  for  adsorbates  have  also 

iv  14 

been  done  for  hydrogen  on  wolfram,  J  hydrogen  on  slllclum,  and  for  the 

Internal  mode  vibration  of  CO  on  nickel.1*  Theoretical  arguments  as  well  as 

detailed  calculations  have  been  given  to  show  that  these  llnewldth 's  are 

1*1-17 

mainly  due  to  pure  dephaslng. 

All  the  above  Indicated  experiments  and  theoretical  considerations  are 
about  the  excitation/relaxation  of  the  vibrational  states  of  a  chemical  bond, 
with  typically  a  frequency  of  2000  cm  1 .  Experiments  of  relaxation  of 
moleoules  embedded  in  the  solid  state  suggest  that  the  line  broadening  is 
mainly  due  to  pure  dephaslng.  There  is  no  reason  to  believe  this  to  be 
different  for  similar  modes  of  adsorbates.1*  The  calculations  in  this  paper 
will  support  these  conclusions. 

However,  for  physlsorbed  systems,  the  adbond  mode  will  have  a  fundamental 
frequency  in  the  region  50  -  500  cm*1,  considerably  lower  than  that  of  a 
chemical  bond.  For  such  low  frequencies  an  experlmantal  analysis  la  difficult 


to  porfora.  We  are  only  aware  of  the  measurement  of  the  llnewldth  of  the 

18 

adbond  vibration  of  CO  on  nickel.  Theoretical  analysis  has  aalnly  been 

19-21 

restricted  to  energy  relaxation.  The  question  we  want  to  answer  in  this 

paper  is  whether  for  these  low-frequency  nodes  pure  dephaslng  still  is  the 

aost  Important  relaxation  process.  It  will  be  shown  that  under  certain 

conditions,  dependent  on  both  the  adbond  and  on  the  details  of  the  reservoir 

spectrua,  Tj-  and  T£-processes  can  become  of  equal  importance,  and  energy 

relaxation  can  even  become  the  more  Important  one. 

19-22 

Following  earlier  work,  In  the  next  section  we  shall  present  a 

theory  of  vibrational  relaxation  of  an  adbond,  with  the  lattice  vibrations  of 
the  substrate  as  a  reservoir.  However,  we  shall  retain  the  terms  describing 
pure  dephasing,  and  we  will  derive  explicit  expressions  for  and  T*.  In 
Secs.  4  and  5,  we  shall  evaluate  these  expressions  numerically  as  a  function 
of  the  fundamental  frequency  of  the  adbond  vibrational  mode.  Finally,  In  Sec. 
6  conclusions  will  be  presented. 


2.  RELAXATION  CONSTANTS 

Consider  a  system  of  an  atom  or  molecule  physlsorbed  on  a  crystal.  It  Is 
assumed  that  that  the  bond  between  the  adsorbate  and  substrate  can  be 
described  as  a  one-dimensional  anharmonlc  oscillator;  only  vibrations 
perpendicular  to  the  surface  are  considered.  Excitation  of  the  adbond  can  be 
accomplished  by  laser  light  In  resonance  with  one  pair  of  vibrational 
states.  Relaxation  occurs  through  the  dynamic  coupling  of  the  adbond  with  the 
lattice  vibrations  of  the  substrate.  The  latter  is  assumed  to  act  as  a 
thermal  bath  or  reservoir.  J 

The  potential  energy  between  the  adatom  and  the  substrate  Is  given  by  a 

one  dimensional  potential  between  the  instanteneous  position  of  the  adatom,  z, 

and  of  the  nearest  surface  atom,  u  (Fig.  1), 

z 

V  -  V(z-z-u)  .  (2.1) 

o  z 

The  origin  is  chosen  as  the  average  position  of  the  surface  atom.  The  average 

of  the  instanteneous  displacement  u  is  therefore  zero.  It  is  assumed  u  is 

z  z 

entirely  determined  by  the  lattice  vibrations  of  the  substrate  (reservoir 
oondition).  Then  the  vibrations  of  the  adatom  are  described  by  the 
Hamiltonian 


H.  •  •  H  ~2  *  <v<*-vV> 

oz 


9 


(2.2) 


where  a  la  the  mass  of  the  adatom  and  <...>  denotes  the  average  over  the 
substrate  lattice  vibrations.  The  Interaction  between  the  adbond  and  the 


lattice  vibrations  Is  given  by 

H1  -  V(z-z0-uz)  -  <V(z-zo-uz)>  (2.3) 

Furthermore,  we  denote  by  Hp  and  p°,  respectively,  the  Hamiltonian  and  thermal 
equilibrium  density  operator  of  the  harmonic  lattice  vibrations  of  the 
substrate.  The  eigenfunctions  and  eigenvalues  of  the  adbond  Hamiltonian  are 
defined  by 

H  |n>  -  Mwln>  .  (2.4) 

a*  n* 

The  complete  system  has  the  Hamiltonian 

H  -  H  ♦  H1  -  H  ♦  H  ♦  H1  .  (2.5) 

o  a  p 

The  state  of  this  system  is  described  by  a  density  operator  p,  which  obeys  the 
equation  of  motion, 

j  CH.p]  .  (2.6) 

Then  the  adbond  can  be  described  by  Its  reduced  density  operator 


o(t)  •  Tr  (p(t)J  , 


(2.7) 


-  8  - 


where  Trp  denotes  the  trace  over  the  phonon  variables. 

19*22 

Following  earlier  work,  an  approxlnate  expression  for  the  reduced 
density  operator  of  the  adbond  can  be  derived  as, 


de(t) 

dt 


1H 


o(t)]  - 

n 


Tr. 


/ 

0 


dt' 


CH1. 


CHI(-f), 


Ppo(t)]3 


(2.8) 


with 


fiX(-t')  -  exp  (-  i  HQt ' )  H1  exp  (i  M0t')  . 


(2.9) 


the  first  term  of  Eq.  (2.8)  defines  the  free  evolution  of  the  adbond  trtille, 
the  second  tern  gives  the  relaxation.  A  further  approximation  can  be  made  by 
neglecting  terms  on  the  right-hand  side  of  Eq.  (2.8)  which  oscillate  with  a 
frequency  different  from  the  free  evolution  of  o(t).  Within  this 
approximation,  usually  denoted  as  the  random-phase  approximation  (RPA),  we 
obtain  for  the  time  derivative  of  the  matrix  elements  of  o(t)  between 
eigenstates  of  H#, 


d«n(t) 

dt 


!  r2  •k<t) 


(2.10) 


and 


d«nm(t) 

dt 


-l  «  o  (t)  -  r“  o  (t) 

bi  rui  ns  no 


(2.11) 


-  9  - 


where  0 _ (t)  -  <nlo(t)|m>,  o„  -  a _ and  *  -  u  -  Let  us  define 

no  n  nn  ruo  n  i 

»■*  -  ^  l  dt’  «">>  "kn<-*,)>  * 

A„  *  1  3  '  «' 

k*n  "  0 


(2.12a) 


(2.12b) 


F 

nm 


/  dt’  «hJ  (0)hJ  ("t* )> 

M2  0  nn  nn 


<HI 


(f)U«)> 

nn 


<iL  (o  )«L  ( -t  • »  -  <hL(  -t  •  )hL  (o )»  , 

no  nn  bus  nn 


(2.12c) 


where  <Q>  •  T  (p°  Q)  is  the  average  of  Q  over  the  phonon  reservoir.  Then  the 
rp  p 

relaxation  constants  are  given  by 


nn  '  ~*kn  for  "  •k 

(2.13a) 

2  '  <*„  *  *„*> 

(2.13b) 

2  *  *•«*„  *  *i  ‘  • 

(2.130 

In  the  RPA,  pure  dephasing  is  given  by  the  diagonal  elements  of  the 

interaction  Hamiltonian.  From  Eqs.  (2.10)-(2. 13)  it  follows  that  indeed  they 

do  not  contribute  to  the  time-evolution  of  the  populations.  They  do  influence 

the  coherences,  with  a  rate  constant  F  and  the  added  the  energy  relaxation 

nm 


rate  constants  A  .  A  oondltlon  for  pure  dephasing  to  occur  is  that  F  is 
n  no 

nonzero.  As  can  be  deduced  tr on  Eq.  (2.12c),  this  means  that 


H1  •  H*  for  n  - 
nn  mb 


(2.14) 


For  example,  It  can  be  shown  that  Eq.  (2.14)  is  not  fulfilled  when  the 

complete  system  is  described  by  harmonic  oscillators.  An  expression  similar 

24 

to  Eq.  (2.12c)  was  given  by  Burns  et  al. 

In  the  rest  of  the  paper  we  shall  use  the  Morse  potential  for  the 
interaction  potential  V(z-z  -u  ), 

O  Z 


V  (z)  -  D{exp(  -2az)  -  2  exp(-az)} 

m 


(2.15) 


The  reasons  are  that  a  Morse  potential  gives  a  fair  description  of  the  van  der 

Waals  bond,  and  that  it  is  to  a  large  extent  amenable  to  analytic  evaluation 

of  the  expressions.  The  lattice  vibrations  are  described  as  a  set  of  harmonic 

19 

oscillators.  Then,  following  the  work  of  Efrima  et  al.,  the  correlation 
functions  of  the  interaction  Hamiltonian  can  be  evaluated: 


<HL(0)Bkl("t,)>  “  °2  e’2“2<u2>  e  l“klt 


b(2)  b(2)  (exp(4o2<u(0)u(-t,)>)-1) 
nm  ki 


♦  4  B<2)  B*!*  {exp(o2<u(0)u(-t')>)-1}  ,  (2.16) 

ns  ki 


-  2  B*2)  B^^  (exp(2a2<u(0)u(-t,)>)-1) 
nn  ki 


■  2  Bnm)  {**P<2«2<u(0)u(-f)>)-1}  , 

where  the  displacement  autocorrelation  function  is  purely  determined  by  the 

(J)  *Jo(z*z' ) 

reservoir.  BVJ'  is  the  matrix  element  of  e  between  the  eigenstates 

nm 

of  H#  (see  the  Appendix). 

3.  MULTIPHONON  CONTRIBUTIONS 

The  Markov  approximation  implies  a  separation  of  the  system  in  two  parts, 
the  reservoir  and  the  subsystem  of  interest  (adbond).  Changes  in  the  state  of 
the  complete  system  then  have  the  Interpretation  of  something  changing  in  the 
adbond  accompanied  by  a  change  in  the  reservoir.  With  in  the  RPA  the  time 
evolution  of  the  reduced  density  operator  of  the  adbond  has  a  very  simple 
Interpretation  (Eqs.  (2.10)  and  (2.11)).  Any  change  of  o(t)  is  accompanied  by 
the  creation  or  destruction  of  one  or  more  phonons  of  the  substrate. 

For  the  Morse  potential  and  a  harmonic  lattice,  the  n-phonon  contribution 
to  the  relaxation  constant  is  given  by  the  combination  of  all  terms  involving 
the  n-th  power  of  the  autocorrelation  function  of  the  surface  atom 
displacement.  Each  of  these  terms  is  singled  out  easily  by  expanding  the 
exponential  functions  in  Eq.  (2.16).  For  a  Debye  spectrum  for  the  phonons, 

the  importance  of  multlphonon  contributions  was  analyzed  by  Jedrzejek  et 

25  26  27 

al.t  *  and  up  to  n  •  2  by  Gortel  et  al.  These  authors  considered  only  the 

master  equation  (Eqs.  (2.10)  and  (2.15)),  and  therefore  analyzed  only  the 

contributions  of  the  multlphonon  contributions  to  the  transition  rate 

constants  a^  (Eq.  (2.12a)).  It  was  found  that  for  transitions  between 

vibrational  states  Involving  energies  less  than  the  Debye  frequency,  the 


s Ingle ’phonon  terns  dominate.  For  transitions  of  energy  higher  than  this 

Debye  frequency,  multi phonon  terms  are  also  important. 

Especially  for  the  motion  perpendicular  to  the  surface,  l.e. ,  u  the 

z 

Debye  spectrum  is  a  poor  presentation  of  the  density  of  states  (DOS)  of  the 
28 

phonon  spectrum.  Recently,  the  transition  rate  constants  for  the 

vibrational  adbond  of  CO  physisorbed  on  Ni  were  calculated  by  Volokitin  et 
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al.  using  the  more  realistic  surface  phonon  spectrum  due  to  Black.  The 

values  they  obtained  were  much  lower  and  in  better  agreement  with  experiment 

than  these  calculated  with  the  Debye  spectrum.  Furthermore,  it  was  found  that 

multiphonon  terms  of  high  order  still  give  a  significant  contribution  to  the 

rate  constants. 

Energy  relaxation  involves  a  change  in  energy  of  the  adbond.  The 
transition  from  level  n  to  level  m  is  accompanied  by  the  creation/anhllatlon 
of  a  phonon,  iftiich  is  depicted  in  Fig.  2.  Multiphonon  transitions  can  be 
represented  by  adding  one  or  more  virtual  adbond  levels.  Pure  dephaslng  (Eq. 
(12c))  occurs  due  to  processes  in  which  the  populations  of  the  adbond  do  not 
change.  These  can  be  interpreted  as  virtual  excitations  of  the  adbond,  iftlch 
then  returns  to  its  original  state.  Obviously  the  minimum  number  of  phonons 
necessary  for  such  a  virtual  excitation  is  two  (Fig.  3).  Therefore,  pure 
dephaslng  has  only  contributions  from  two-phonon  and  higher  multiphonon 
processes. 

From  the  notion  that  pure  dephasing  is  at  least  a  two-phonon  process,  and 
from  the  previous  discussion,  one  expects  that  pure  dephasing  will  become  an 
important  relaxation  process  whenever  the  multlphonon  contributions  to  the 
transition  rate  constants  are  important.  The  analysis  in  the  rest  of  the 
paper  will  support  this  hyphothesls  to  a  large  extent,  although  it  will  also 


be  found  that  the  Intrinsic  differences  between  pure  dephasing  and  energy 
relaxation  do  influence  their  relative  Importance. 

Let  us  denote  the  phonon  density  of  states  by  p(w),  with 

m 

f  d <d  p(u)  •  1  (3.1 ) 

0 

Then,  the  displacement  autocorrelation  function  Is  given  by, 

<u(0)u(-t')>  •  ~  /  d«  {WiiJMje’1"1'  ♦  ^coJe1"1'}  (3.2) 

2M  o  “ 

with 

n(w)  -  (exp( hu/kgT)  -  1)  1  ,  (3-3) 


which  gives  the  number  of  phonons  with  frequency  w  as  a  function  of  the 
temperature  T;  kg  is  the  Boltzmann  constant. 

With  these  definitions  the  two-phonon  contribution  to  the  pure  dephaslng 
rate  constant  becomes  for  a  Morse  potential  (Eq.  (2.12c)), 


Re(F  ) 

nm 


^2aW> 


.m-n.2  Hi 
k  M2 


2w  - 

(SlU)2 

l2M  ' 


«  /  dw  (&£y~)2  n(w)[n(u)  *  1]  • 

0  w 


(3.H) 


I 


-  in  - 


where  D,  a  and  k  are  defined  in  the  Appendix.  The  inag inary  part  of  the  rate 

constants  deteralnes  a  level  shift  and  will  not  be  considered  here. 

The  energy  relaxation  contribution  to  the  absorption  llnewidth  of  a 

transition  n  ->  a  Is  given  by  Re(A„  ♦  A*)  -  Re(A«  *  A  )  (Eq.  (2.13c)).  Re(A  ) 

nans  n 

can  be  written  as  (Eq.  (2.12)) 
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For  the  latter  rate  constants  the  1  -phonon  contribution  is 
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The  two-phonon  energy  relaxation  is  given  by. 
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The  four  terms  in  Eq.  (3.7)  correspond  to  the  four  two-phonon  processes 
denoted  by  II-V  in  Figure  2. 

From  Eqs.  (3.6)  and  (3.7)  we  can  conclude  that  the  energy  relaxation 
constants  strongly  depend  on  the  details  of  the  phonon  spectrum  used.  In 
particular,  the  single-phonon  term  Is  directly  proportional  to  the  DOS  at  the 
transition  frequency.  The  pure  dephasing  constant  Is  less  dependent  on  the 
phonon  spectrum  because  It  involves  an  Integral  over  the  square  of  p(u).  In 
the  next  two  sections  we  shall  evaluate  the  expressions  Eqs.  (3.*0.  (3*6)  and 
(3.7)  for,  respectively,  a  Debye  spectrum  and  a  DOS  representing  a  surface 
phonon  spectrum. 


*.  DEBYE  SPECTRUM 

For  a  Debye  spectrum  we  have 
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Let  us  define  the  integral 


1  x2ex 
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0  (e  -1 ) 


Then  the  two-phonon  contribution  to  the  pure  dephaalng  for  a  Debye  phonon 
spectrun  can  be  written  as, 
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It  is  easy  to  show  that  ID(y)  •  constant  -3.3  for  (y  ♦  -).  Then  we  obtain 

Re(F  )  -  0(T3)  for  T  •*  0. 
nn 

Substitution  of  Eq.  (M.1)  in  Eq.  (3.6)  gives  for  the  single-phonon  energy 
relaxation 


i(1) 

no 


D2  .***<«*>  (B*'*)2  (n-.)2  (1 

nn 


6* 


n*in*1 .  .a  . 
2k  1  'M  ; 


(M.h) 


*«  c'  * 
D 


.(1) 


for  «  i  «D.  For  m  >  we  have  a^  -  0.  The  two-phonon  contribution  is 
(Eq.  (3.7)), 


a(2)  - 

nn 


Mu>. 


■>’  .-**«*>  {p  ^  .0  , 

""D  B  8 


(*.5) 


-  17  - 
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The  expression  for  the  integral  E  depends  on  the  values  of  «D  and  uf 
for  |ej  >  2^  we  have 
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with  sinilar  expressions  for  <  0. 


For  T  •*  0,  x..^  converges  to  a  lloiting  value.  This  aeans  that 

V  u  uB 


at  low  teaperature  all  n-phonon  contributions  to  the  energy  relaxation  remain 


finite.  Therefore,  in  the  llBlt  T  ♦  0,  energy  relaxation  will  always  dominate 


pure  dephasing. 
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We  ere  interested  In  en  estimate  of  the  different  contributions  to  the 

total  relaxation  rate  Eq.  (2.13c).  In  order  to  do  this  we  shall  use  values 

for  the  parameters  which  represent  realistic  systems  like  Ar/KCl  and  C0/N1. 

That  Is,  we  use  a  number  of  adbond  levels  k  •  25.  Different  systems,  i.e. , 

with  a  different  fundamental  frequency,  are  then  represented  by  choosing 

different  values  for  both  a  and  D,  such  that  k  »  25.  The  values  for  v  and  M 

26 

are  taken  in  accordance  with  the  system  CO/Cu.  In  the  low “temperature  limit, 
Xp  <  1 ,  the  energy  relaxation  rate  constant  Re(A^Ao)  for  the  1  ♦  0  transition 
can  be  approximated  by  ^a1Q.  Also,  for  the  1  ♦  0  transition  (n,m  <<  k),  the 
n,s-dependent  factors  have  a  simple  approximate  expression.  With  these 
substitutions  and  approximations,  the  relaxation  oonstants  become 


It  Is  noted  that  for  «10  >  Wp.  aj^ 


0,  but 


-  0,  and  that  for  «10  >  2^,  a^'  • 
that  the  pure  dephaslng  constant  is  only  weakly  dependent  on  the  fundamental 
frequency. 

From  Eq.  (4.7)  the  following  estimates  of  the  relative  importance  of  the 
various  contributions  at  low  temperature  can  be  made  for  0  <  u10  <  «D: 


(4.8a) 


(4.8b) 


For  <  2*>p,  the  single-phonon  energy  relaxation  is  identically  zero, 

but  Eq.  (4.8b)  remains  valid. 

From  Eq.  (4.8),  it  follows  that,  for  a  Debye  spectrum,  the  pure  dephaslng 
rate  constant  is  such  smaller  than  the  single-phonon  energy  relaxation,  and  of 
about  the  same  magnitude  as  the  two-phonon  energy  relaxation.  The  two-phonon 
contribution  to  the  energy  relaxation  is  much  smaller  than  the  single-phonon 
term,  and  the  three-phonon  term  will  again  be  an  order  of  magnitude  smaller 
than  the  two-phonon  term.  Therefore,  when  compared  to  three-phonon  or  higher- 
order  terms  in  the  energy  relaxation  rate  constant,  the  pure  dephaslng  will  be 
the  most  important  relaxation  process. 


5.  SURFACE  PHONON  SPECTRUM 


In  the  aodel  used  in  this  paper,  the  coupling  between  the  adbond  and  the 

phonon  reservoir  is  nedlated  by  the  notion  of  the  surface  aton  perpendicular 

to  the  surface.  Especially  for  this  direction  of  the  notion,  the  theoretical 

28 

calculations  of  Black  for  a  Ni (111)  surface  show  a  large  difference  in  the 
spectrun  of  the  surface  phonons  conpared  to  that  of  the  bulk  phonons.  The 
nost  prominent  characteristic  of  this  spectrun  is  that  the  DOS  consists  of  a 
single  peak  at  frequency  eQ  •  24.2  THz.  Besides  this  peak,  there  is  sons 
extra  density  in  the  tails  between  frequencies  of  60  cn  and  300  cn  (as 
visible  in  the  picture  given  by  Black).  We  want  to  evaluate  the  expressions 
for  the  relaxation  constants  for  such  a  phonon  density  of  states.  In  order  to 
do  this,  we  approximate  the  calculated  spectrun  by  the  slnple  analytical  form 


p(w) 


— y-  u*  exp(-|e 1  -w|ln2/A) 


2"oA 


(5.1) 


1  2 

where  A  -  1.45  THz  -  The  u  -dependence  is  chosen  in  order  to  obtain 

the  correct  Uniting  behavior  for  »  ♦  0;  the  exponential  gives  a  reasonable 
representation  of  the  spectrun  calculated  by  Black. 

Again,  we  consider  the  1  ♦  0  transition  in  the  low-teaperature  Unit. 
Substitution  of  the  DOS  (Eq.  (5.1))  into  Eqs.  (3*4),  (3.6)  and  (3.7)  gives 
the  required  rate  constants.  After  rearrangement  they  becose 
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where  z  •  «/«o,  and  Ig(-^)  is  given  by 
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±  (f  "  x  *  ?  *^>)  «*P(“  ln2  |2-x|  *0/A) 

for  x  i  1.  The  plus  sign  applies  for  x  2  2  and  the  alnus  sign  for  11x12 

(x  *  •>/»»).  It  follows  that  the  energy  relaxation  rate  constants  are,  In 
na  o 

this  low-temperature  approx lmat Ion,  independent  of  the  temperature.  Pure 
dephaslng  is  strongly  temperature  dependent,  through  the  exponential 
functions. 

The  Norse  potential  parameters  representing  soae  real  physisorbed  systems 

26 

are  such  that  the  number  of  levels,  k,  is  more  or  less  constant.  To  compare 
the  energy  relaxation  and  pure  dephaslng,  we  plot  in  Fig.  4  the  relaxation 
constants  as  a  function  of  the  fundamental  frequency,  w1Q,  of  the  vibrational 
adbond.  The  constants  defining  the  adbond  (k,  v ,  H)  are  the  same  as  those 
used  In  Sec.  4. 

All  quantities  are  given  in  units  of  the  characteristic  frequency  u  of 

o 

the  phonon  spectrum.  Ve  see  that,  as  expected,  the  n-phonon  energy  relaxation 

rate  constant  has  a  maximum  at  frequency  nw  .  For  •  u  it  decreases 

o  10  o 

exponentially  with  Im^-w  |.  The  pure  dephaslng  constant  is  much  less 

dependent  on  the  transition  frequency.  The  largest  effect  is  due  to  the 
d2V(*-z  ) 

increase  of - = -  |  with  increasing  *10,  and  not  to  a  resonance  effect 

dz  z-lo 

as  for  the  energy  relaxation. 

It  is  to  be  expected  that  the  maxima  in  aj^  will  gradually  decrease  with 
increasing  n.  Then,  for  transition  frequencies  larger  than  n«Q  with  n  on  the 
order  3  or  U,  pure  dephaslng  will  become  a  more  important  process  than  energy 
relaxation.  A  notable  difference  with  the  Debye  spectrum  is  that  the  single- 
phonon  and  two-phonon  contributions  to  the  energy  relaxation  are  still  larger 
than  the  pure  dephaslng  even  if  >  nuQ  (n  -  1,2,  respectively).  For  the 


Debye  spec  true  the  n-phonon  tern  Is  exactly  zero  for  >  nup.  This  is  due 

to  p(«)  now  being  non-zero  for  u  >  mQ.  Therfore,  the  precise  behaviour  of 

*10^  for  *10  >  ”“0  dePend8  strongly  on  the  exact  form  of  p(w).  Since  the 

28 

difference  between  the  calculated  spectrum  and  a  Lorentzlan  peak  (used  here) 
or,  for  example,  a  Gaussian  function  cannot  be  seen  for  t>  2  2«qP  the  values  of 
the  n-phonon  relaxation  constants  beyond  ■  (n*1)«0  are  not  reliable,  and 
no  conclusion  should  be  drawn  regarding  their  behavior  at  those  frequencies. 

6.  DISCUSSION 

It  is  well  known  that  pure  dephasing  is  a  relaxation  process  which  is  due 

to  anharmonicltles  In  the  system.  In  the  model  used  In  this  paper, 

anharmonlclty  is  Introduced  in  the  active  mode  (vibrational  adbond),  and  its 

magnitude  Is  fixed  by  the  choice  of  a  Horse  potential  to  represent  this  bond. 

Another  approach  can  be  to  take  Into  account  the  anharmonlclty  in  lattice 

vibrations.^  We  assume  that  this  can  be  neglected  with  respect  to  the  effect 

of  the  adbond  anharmonlclty.  Therefore,  it  should  be  kept  in  mind  that  the 

17  11 

calculations  presented  here  and  elsewhere  give  values  for  the  pure 
dephaslng  constants  which  can  be  strongly  model  dependent.  Energy  relaxation 
Is  already  present  in  a  harmonic  system  and  does  not  change  dramatically  upon 
Introducing  anharmonlclty. 

It  is  found  that  for  low  temperature  the  pure  dephaslng  constant  goes  to 
zero  as  (T  ♦  0)  (Eqs.  (A.7)  and  (5.*);  see  also  Ref.  17).  Apparently  this 
result  does  not  depend  on  the  details  of  the  phonon  spectrum.  It  does  depend 
on  the  kind  of  relaxation  process  studied,  l.e.,  on  the  fora  of  the  coupling 
between  the  active  mode  and  lattice  vibrations.  This  can  be  compared  to  the 
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results  of  Skinner  and  Hsu  who  find  a  T -dependence  for  the  relaxation  of  the 
electronic  excitation  of  a  molecule  embedded  in  a  crystal.^1 

The  calculations  presented  in  Secs.  3“5  clearly  show  that  resonance  is  an 
important  effect  in  determining  the  magnitude  of  the  energy  relaxation  rate 
constants  (Eqs.  (3.6)  and  (3.7)).  The  n-phonon  contribution  is  large  whenever 
all  phonons  involved  are  of  i  frequency  at  which  the  DOS  is  large.  For  the 
Debye  spectrum  this  means  that  the  n-phonon  term  has  a  maximum  for  the 
transition  frequency  smaller  but  of  the  order  nw^.  For  a>1Q  >  nu^,  that  term 
is  identically  zero.  For  the  surface  phonon  spectrum,  each  n-phonon  term  in 
the  rate  constant  has  a  maximum  at  nu0  and  decreases  strongly  for  |wj0-n<i»o|  * 
0.  The  results  also  show  that  on  going  to  a  higher  multiphonon  process,  the 
maximum  value  of  the  relaxation  constant  decreases  with  Increasing  n  (see  also 
Refs.  25-27). 

The  lowest-order  contribution  to  the  pure  dephasing  is  a  two-phonon 
process.  Because  it  Involves  only  a  single  adbond  level,  the  resonance  effect 
is  now  that  the  largest  contribution  to  the  relaxation  constant  cooes  from 
virtual  excitations  with  frequency  uQ.  Since  there  are  no  other  restrictions 
on  the  position  of  the  virtual  level,  the  relaxation  constant  will  hardly 
depend  on  the  fundamental  frequency.  Its  magnitude  should  therefore  be  of  the 
order  of  the  two-phonon  energy  relaxation  term  at  its  maximum.  On  the  other 
hand,  it  will  roughly  be  lower  by  a  factor  of  n(«0>  than  this  term,  because  it 
also  Involves  the  absorption  of  a  phonon.  It  Is  this  factor  which  introduces 
the  strong  temperature  dependence  of  the  pure  dephaslng  (Fig.  4). 

This  analysis  shows  that  at  high  temperatures,  where  ntw^)  or  n(w0)  >  1, 
the  two-phonon  pure  dephaslng  term  will  be  of  the  order  of  the  maximum  value 
of  the  two-phonon  energy  relaxation  term.  For  lower  temperature  (n(«o),  n(u>D) 


<  1),  the  two-phonon  pure  dephasing  will  be  somewhat  smaller.  Then,  one 

expects  pure  dephasing  to  be  the  most  Important  relaxation  process  when  the 

transition  frequency  of  the  active  mode  is  much  larger  than  the  characteristic 

frequency  («o  or  w^)  of  the  phonon  reservoir.  This  is  in  accordance  with  the 
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result  of  Gadzuk  and  Luntz  and  Hutchinson  and  George,  who  considered  the 
relaxation  of  an  Internal  vibrational  mode  of  an  adsorbed  molecule.  Although 
the  coupling  mechanism  is  different  there,  the  conclusions  are  also  valid  for 
the  relaxation  of  an  electronic  excitation.^1  However,  for  the  low-frequency 
adbond  vibrations  of  a  physisorbed  molecule,  the  energy  relaxation  is  mainly 
due  to  single-phonon  and/or  two-phonon  contributions.  Then  energy  relaxation 
will  give  the  largest  contribution. 
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APPENDIX:  Matrix  Element 
-  - - — nm 

The  Morse  potential  is  given  by 

V*-zo)  -  D{exp(-2a(z-z0>)  -  2  exp(  -o(x-*0)))  . 

Averaging  V(z-zo-u)  over  the  lattice  vibrations  (i.e. , 
again  a  Morse  potential, 

<V(z-zft-u)>  -  V'(z-z')  , 

no  no 

with  renormalized  parameters 


D'  -  D  exp(-a2<u2>) 

o'  •  a 


z '  *  z  ♦  |a<u2> 
O  0  2 


The  adbond  Hamiltonian  is  defined  with  the  renormalized 
potential  energy.  Let  us  its  decrete  eigenfunctions  by 
corresponding  eigenvalues  are 


w2  2  ,  - 
n  2p  2} 


where 


(A  1 ) 

over  u)  results  in 

(A2) 

(A3) 

Morse  potential  as 
|n>.  Then  the 

(AH) 


fcnv.-v„-> .TJ.W v-^v.  a w v- kj, v  v-v.  ■> ■> 


I 


The  Integer  part  of  k  represents  the  number  of  bound  states.  From  Eq.  (AH)  one 


directly  obtains  the  fundamental  frequency 


“10  ■  *r- (k‘”  • 


The  following  matrix  elements  can  be  calculated  as 


19,32 


•flf  I  «7  •  ^ 

B(1)  -  <n\m  0  l»>  -  J.  f(2k-2n-1)(2k-2m-1).  1/2  r(n.Ur(9u^ 

Bnm  <nl*  '  2k  {n!m!r(2k-n)r(2k-m) }  r(  1)r(2k'®)  (A7) 


B*2)  -  <n|exp(-2a(z-z')|m>  -  bJJ*  ( I4 j* )  f  (a8) 


both  for  n  >  (m-1). 


We  need  the  matrix  elements  of  the  Interaction  Hamiltonian  (Eq.  (2.12)). 


Substitution  of  the  explicit  form  of  H  (Eq.  (2.3))  shows  that  this  involves 


itrlx  elements  of  expC-Jo(z-zQ)] ,  for  J  •  1,2.  That  Is,  we  need  matrix 


elements  of  both  terms  of  the  Morse  potential  (Eq.  (A1))  between  the 


eigenstates  of  the  renormalized  Morse  potential  (Eqs.  (A2)  and  (A3)).  We  note 


that  exp[-Ja(z-z  )]  -  expC-Ja(z-z')]exp[-Ja(z'-z  )] ,  and  that  (z’-z  )  -  |a<u2> 

O  O  O  O  0  0  2 


Is  constant.  Then  from  Eqs.  (A7)  and  (A8)  we  immediately  obtain  the  result 


B**J  -  <njexp(-ja(z-z  )|m>  -  e 


-  32  - 


with  *  •  1,2. 

This  result  was  used  in  the  evaluation  of  Eq.  (2.16).  We  stress  that  the 
correction  so  obtained  (Eq.  (A9)  versus  Eqs.  (A7)  and  (A8))  is  of  the  same 
order  of  aagnitude  as  the  corrections  obtained  directly  froo  using  the 
renormalized  Morse  potential  Instead  of  VB(z-z0)  itself  in  H#. 
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FIGURE  CAPTIOUS 

Figure  1.  Schematic  drawing  of  an  adsorbed  atom  at  position  z  a bore  a  surface. 
The  nearest  surface  atoa  has  a  time -dependent  displacement  u 
perpendicular  to  the  surface. 

Figure  2.  The  single-phonon  and  two-phonon  contributions  to  the  energy 

relaxation  for  the  transition  from  level  n  to  level  ■,  with  m  >  <■>  . 

n  b 

Each  arrow  denotes  the  absorption  (upwards)  or  eelssion  (downwards) 
of  a  phonon.  The  energy  of  each  phonon  follows  flrom  the  length  of 
the  arrow.  The  dotted  lines  indicate  virtual  adbond  levels.  I  is 
the  slr^le-phonon  process,  and  there  are  four  two-phonon  processes 
indicated  by  1I-IV.  For  the  downward  transition  n  ♦  m,  process  III 
Is  not  possible,  and  for  the  reverse  upward  transition,  process  II 
does  not  contribute. 

Figure  3»  Miltiphonon  contributions  to  pure  dephaslng.  Depicted  are  both  two- 
phonon  processes  and  one  three -phonon  process.  For  the  meaning  of 
the  symbols  see  Fig.  2.  Note  that  pure  dephaslng  terms  Involve  only 
one  (real)  level  of  the  adbond,  and  that  single-phonon  proceses  are 
impossible. 

Flgurm  4.  The  logarithm  of  the  various  n-phonon  contributions  to  the 

relaxation  constants  versus  the  fundamental  frequency  of  the 

vibration.  The  low -temperature  approximation  is  used  (Eqs.  (5.2)- 

(5.4)).  All  quantities  are  given  in  units  of  the  characteristic 

(2) 

frequency  of  the  surface  phonon  spectrum.  The  pure  dephaslng,  F^0  , 


v*  •/  s' 


-V.y..- 


is  plotted  for  two  values  of  the  temperature.  Note  that  the  low- 

**10 

temperature  approximation  does  not  make  sense  for  — <  r~~. 
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